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The first VLBI images of H2O maser emission in the Circinus Galaxy AGN show both an 
accretion disk and outflow 0.1 to 1 pc from the central engine. The maser traces a warped, 
edge-on accretion disk between radii of about 0.1 and 0.4 pc that is bound by a 1.3 x 10 6 
M© central mass. The rotation curve is somewhat flatter than a Keplerian rotation law and 
is consistent with a disk mass on the order of 10 5 M©. Clumpy substructure may display a 
kinematic signature consistent with spiral arms. Away from the disk, a second population of 
water masers traces broadly a bipolar, wide-angle outflow that contains (bullet-like) clumps 
ejected from a region centered on the central engine and < 0.1 pc in radius. Out to a radius 
of ~ 0.4 pc, the warp of the accretion disk appears to channel the outflow, the orientation of 
which on the sky coincides with the orientation of the known kiloparsec-scale ionization cone. 
Beyond this radius, the flow crosses the disk and truncates it. This suggests that the current 
accretion event and associated disk-outflow geometry have a lifetime on the order of 10 7 years. 



1. Introduction 

The Circinus galaxy is one of the nearest Seyfert II galaxies that hosts an extragalactic 
H2O maser. (Freeman et al. 1977 estimate a distance of 4 Mpc.) The ~ 10 42 erg s _1 (2- 
10 keV) central engine is obscured at energies below 10 keV by a large gas column, nn ~ 
4 x 10 24 cm~ 2 (Matt et al. 199£). The X-ray spectrum also exhibits a prominent Fe Ka 
line that probably comprises reflected radiation (Matt et al. 1996). These characteristics 



suggest some similarity to another maser-host galax y, NGC 1068. 

The Circinus AGN exhibits an ionization cone ( Marconi et al. 1994 ) within which 
Veilleux & Bland-Hawthorn (1997) observed linear optical filaments and compact knots 
reminiscent of Hcrbig-Haro objects. The cone opening angle is at least 90°, at a mean 
position angle (PA) of roughly 2 90°. The AGN also d rives a kpc-scale nuclear outflow 
betrayed by bipolar radio lobes ( Elmouttic et al. 1998) at PA ~ 29 5°, which are largely 
align ed with the minor a xis the galactic HI disk ( Jones et al. 1999 ) and a nuclear 12 CO 
ring ( Curran et al. 199S ). 

There is strong evidence that at least several extragalactic H2O masers trace edge-on 
accretion disks bound by central engines > 10 6 M Q . These masers display (1) roughly 
linear structure in projection on the sky, (2) emission symmetrically bracketting the 
galactic systemic velocities, and (3) declining rotation curves that indicate more or less 
Keplerian differential rotation. Discovery of similar symmetry in the Circinus H2O maser 
spectrum (Nakai et al. 1995; Greenhill et al. 1997) led to speculation that there too the 
masers trace an accretion disk (Greenhill et al. 1997). Though the observation of an 
accretion disk in Circinus confirms that prediction, the discovery of qualitatively new, 
nondisk emission is equally interesting. In NGC 1068, disk emission is accompanied by 
secondary maser emission ^> 1 pc downstream in a jet, where material in the narrow line 



region deflects it (Gallimore et al. 1996). However, the nondisk emission in Circinus lies 
< 1 pc from the central engine and is more closely related to it. 



2. Observations and Data 

We observed the Circinus H2O maser three times for 18 hours each in 1997 (June & 
July) and 1998 (June) with four stations of the Australia Telescope Long Baseline Array, 
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obtaining a half-power bcamwidth of about 2x4 milliarcseconds (mas) . We recorded two 
16 MHz bandpasses (~ 215 km s _1 ) at each station and correlated the data with the S2 
processor in Marsfield, NSW . Following Fourier inversion of the correlation functions, 
each spectral channel corresponded to ~ 0.21 km s _1 , which provided at least three 
channels across the half-power full- width of each spectral component. 

We calibrated amplitudes, delays, and phases with standard VLBI techniques and the 
AIPS package. When possible we fitted a time-series of total-power spectra for each 
station to a calibrated template spectrum and computed multiplicative gain factors that 
would calibrate visibility amplitudes in units of Janskys. This calibration corrected for 
variations in antenna gain and in atmospheric attenuation and emission, with a < 10% 
relative accuracy (station-to-station, moment-to-moment) and 30% accuracy in absolute 
terms. When during one 1997 track the total-power signal was too weak for template 
fitting, we adjusted the gains to maintain constant peak visibility amplitude in the cross- 
power spectra. We used scans of 0537-448, 1144-379, 1424-418, and 1921-293 to estimate 
(and correct for) instrumental delays and fringe rates at better than the 2 ns and 2 mHz 
level (1998) or 5 ns and 5 mHz level (1997). We also estimated the stations' complex 
bandpass responses from observation 1424-418 and 1921-293. 

We estimated a new astrometric position for the masers by analyzing the time variation 
in fringe rate for the 565 km s _1 line. (We use the radio astronomical definition of Doppler 
shift.) The new position, a 20 oo = 14H3 m 09f95 ± S .02, <5 2 ooo = -65°20'21"2 ± O'.'l, is 
the most accurate estimate so far for the AGN. 

We self-calibrated the emission within a few km s _1 of 565 km s _1 , applied the cal- 
ibration to each spectral channel, and created deconvolved synthesis images. The noise 
in the images (la) was 0.025 - 0.045 Jy, depending on the epoch and channel; the 1998 
June observations were the most sensitive and best calibrated, because of equipment 
upgrades and a (coincidental) 40 Jy flare that made possible particularly accurate am- 
plitude and phase calibrations. We fitted 2-D Gaussian model brightness distributions 
to each emission component in the deconvolved images that was stronger than 6<r. 

At each epoch, the emission near 565 km s _1 was strongest. However, because spectral 
lines come and go from epoch to epoch, we superposed maps of the emission component 
positions from each epoch, so as to trace the underlying dense molecular gas as completely 
as possible with the available data (Figure 1). We registered the maps with 0.2 mas 
accuracy by comparing the measured positions of the lines at 472 and 526 km s _1 . (To 
achieve the stated noise levels, the 1997 data were smoothed using a moving average of 
up to six frequency channels, which "blurs" the map in crowded fields. The 1998 data 
were Hanning smoothed.) 



3. Results 

The maser emission comprises complexes of spectral lines that are redshifted and 
blueshifted with respect to the systemic velocity of the galaxy. In 1998 June, we detected 
emission between 214.8 km s _1 and 676.5 km s _1 , a broader interval than the roughly 
250 to 650 km s" 1 range reported by Nakai et al. (1995) and Greenhill et al. (1997) in 
total-power observations. The mean velocity, 446 km s _ , is consistent with the systemic 



HI velocity of the galaxy, 438 ± 2 km s _1 (Freeman et al. 1977), though emission on 
scales > 100 pc dominates measurement of the latter. 

At each epoch, the sky distribution of maser emission may be divided into three pop- 
ulations (see Figure 1): (1) a thin, gently curved, S-shaped locus of highly redshifted 
and blueshifted emission arcs to the southwest and northeast, respectively (aka "high- 
velocity" emission), (2) emission close to the nominal systemic velocity of the galaxy 
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(aka "low-velocity" emission) that lies between the high-velocity arcs, and (3) emission 
that is modestly Doppler shifted and broadly distributed in knots that lie north and west 
(redshifted) and south and east (blueshifted) of the low-velocity emission. 

3.1. The Warped Disk 

To the position and velocity data for high and low-velocity masers, we have fitted a model 
edge-on disk with smoothly varying position angle as a function of radius (Figures 2 & 
3). The disk outer radius is ~ 0"02 (0.4 pc) and the thickness is < 0''002 thick. The 
orbital speed along the inner edge of the disk, at a radius of ~ 0.1 pc, is 237 km s^ 1 . The 
mass enclosed is 1.3 ± 0.1 x 10 6 M©, assuming circular motion, and the model systemic 
velocity is 451 ± 10 km s . (Uncertainties reflect formal errors.) Among the redshifted 
high-velocity masers, the peak rotation velocity as a function of impact parameter from 
the dynamical center, b, traces a rotation curve that declines as approximately 5-°- 4 ± 05 
(Figure 3), from which we infer a disk mass on the order of 10 5 M© between 0.1 and 0.4 
pc. The disk is edge-on at a radius of 0.1 pc, but if the disk becomes inclined toward the 
outer radius, then the disk mass would be larger. 

The model derives its strongest support from the following observations: (1) the sky 
distribution of high- velocity masers is highly elongated, roughly symmetric, and perpen- 
dicular to the approximate axis of the known ionization cone, (2) the apparent rotation 
curve is nearly Keplerian, (3) the most highly blue-shifted emission, innermost redshifted 
emission, and the low-velocity emission lie along a line on the sky, and (4) a single 
position-velocity gradient extends from the high-velocity blueshifted emission through 
the low-velocity emission to the high-velocity redshifted emission. The latter two obser- 
vations are strong signatures of an edge-on disk in which the low- velocity emission lies 
close to the inner radius of the disk, as it does in NGC4258 ( Miyoshi et al. 1995| ). 



Together with the shallow rotation curve, apparent clumpy substructure suggests that 
self-gravity is important in the disk. Estimates of the Toomre Q-parameter suggest 
that the disk is at best only marginally stable, due in part to its slow rotation speed 
(cf Maoz (1995) in discussion of NGC4258). At the inner radius, where ft is greatest, 
Q ~ OApg 1 hzl, where pg is H2 particle density in units of 10 9 cm -3 and h^i is the 
disk thickness in units of 0.01 pc. (Q = c s VL/nGY>, where c s is the speed of sound, 
~ 1 km s _1 at 400 K, is the angular rotation speed, and £ is the surface density.) 
Furthermore, the kinematic signature of at least two clumps among the redshifted high- 
velocity masers may indicate that spiral arms are present, specifically, the decline in 
line-of-sight velocity with increasing impact parameter. We infer that these clumps are 
elongated along the line of sight and cross the disk midline (i.e., the diameter perpendic- 
ular to the line of sight) at an angle. The angle is roughly the same in each case, which 
is suggestive of ordered structure, though the putative spiral may be fragmented. 

3.2. The Outflow 

A flow from the vicinity of the central engine is traced by: (1) emission between 300 and 
450 km s _1 south and east of the dynamical center of the disk, and (2) emission between 
450 and 580 km s , north and west of the dynamical center. This model demands that 
individual clumps are sufficiently luminous to be observable despite their great distance, 
but this is reasonable. A typical maser in the putative outflow is 0.1-0.5 Jy at a distance 
of 4 Mpc and would appear to be 2 — 8 x 10 4 Jy were it at the distance of the W49N star 



forming region, which is not unusual (Liljestrom et al. 1989). In contrast, it is unlikely 



that the maser emission arises in the accretion disk because of the irregular distribution 
of line-of-sight velocities on the sky, and because were the emission to arise in a disk, it 
would not originate where amplification paths (as dictated by geometry) are optimal. 
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The segregation of red and blueshifted outflowing clumps on the sky suggests that the 
flow is inclined with respect to the line of sight, while the distribution on the sky suggests 
a wide-angle outflow. However, our census of outflowing clumps is incomplete because 
the inter-clump medium is probably ionized, as in narrow-line regions, and readily able 
to obscure emission from the far side of the outflow. (The known emission is close to our 
detection limit and an optical depth of a few requires only readily attainable emission 
measures of 5 x 10 9 pc cm -6 .) Because H 2 maser emission requires H 2 densities ~ 10 9 
cm -3 , the observed clumps are probably high-density "bullets" immersed in a thinner 
medium. Comparable interclump densities would imply a prohibitively large mass loss 
rate (3> 10 M©yr -1 ) and unusually small accretion efficiency. 

In principle the outflowing clumps may be ejected ballistically from the vicinity of 
the central engine (at radii < 0.1 pc) or shorn from the observed molecular disk (e.g., 



Emmering, Blandford, & Shlosman 1992, Kartje, Konigl, & Elitzur 1999). In the absence 
of acceleration or proper motion data, ejection may be the preferred (and simplest) 
model for three reasons. First, some clumps lie close in projection to the rotation axis 
of the accretion disk across a range of impact parameters (i.e., from the disk inner 
radius to at least twice the disk outer radius) and velocities. Second, the line-of-sight 
velocity and position data do not betray rotation. Third, the physical conditions in 
the maser clumps are consistent with adiabatic expansion from conditions estimated for 
broad line clouds (i.e., densities of 10 10 - 10 12 cm -3 and temperatures of 10 4 K, as in 



Brothcrton et al. 1994[) 



The disk structure and rotation curve dissolve at a radius of ~ 0.4 pc, where the wind 
appears to disrupt the accretion disk, whose density presumably declines with increasing 
radius. In the model of Neufeld, Maloney, & Conger (1994) maser excitation depends on 
irradiation of molecular gas by hard X-rays. The proposition that the disk extends no 
farther than 0.4 pc from the central engine is strengthened by the circumstantial evidence 
that the outflow-borne masers all lie outside the shadow cast by the truncated disk. This 
would not be true for a broader disk. However, if the outflow subtends almost 90 % of 
Air str as seen from the central engine, then the accretion disk may be starved such that 
it will be exhausted on the order of 10 7 yr, assuming a mass of 10 5 M , 10% accretion 
efficiency, and 10% emission fraction for hard X-rays. 

The position angles at which the outflow is free of occultation by the 0.4 pc radius 
disk coincide with the limits of the kiloparsec scale ionization cone observed west of 



the nucleus ( Vcillcux fc Bland-Hawthorn 1997 ). The southern edges of the redshiftcd 



maser outflow and the ionization cone, both lie at position angles of ~ —120°. The 
northern edge of the maser outflow, at 16°, coincides roughly with the position angle 
of the northernmost [OIII] filament (~ —20°) or a northeastern blue knot (~ 10°). In 
addition, the mean axis of the maser outflow (—52°), corresponds well to the orientation 
of the domi nant [OIII] filament ( ~ —50°) and a radio hotspot, observed with arcsecond 
resolution ( |Elmouttic et al. 1998 ), which could mark a radio jet, though none has yet 



been observed directly. 

This work has been conducted in collaboration with R. S. Booth' 1 ), S. P. Ellingsen' 2 ) , 
J. R. Herrnstein' 3 ), D. L. Jauncey^, P. M. McCulloch^, J. M. Moran' 5 ), R. P. Norris, 
J. E. Reynolds, and A. K. Tzioumis' 4 ). 

^Onsala Space Observatory; ^ University of Tasmania; ^'Renaissance Technologies; 
(^Australia Telescope National Facility; ( 5 ) Harvard-Smithsonian Center for Astrophysics 
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Figure 1. Sky distribution of H2O maser emission at three epochs (superposed). The synthe- 
sized beam for each epoch is indicated. Symbol shading indicates line-of-sight velocity for each 
spectral channel plotted. Error bars indicate total position uncertainties (la). 
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Figure 2. Position- velocity diagram annotated to highlight the proposed disk and outflow 
components. Dashed rotation curves correspond to the midline of the disk. The steep diagonal 
line corresponds to the near side of the disk, at the inner radius. The outflow component is 
Doppler shifted with respect to the systemic velocity by on the order of ±100 km s _1 . 
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Figure 3. Model of the disk-outflow structure in the inner 1 pc of the AGN overlayed on 
the maser emission distribution. Dashed lines indicate a warped edge-on disk. Outward facing 
arrows indicate a wide-angle outflow. Regions without a direct line of sight to the central engine 
are shadowed by the disk. 
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